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We investigate a two-component atomic Fermi gas with population imbalance in the presence of 
Rashba-type spin-orbit coupling (SOC). As a competition between SOC and population imbalance, 
the finite-temperature phase diagram reveals a large varieties of new features, including the expand- 
ing of the superfluid state regime and the shrinking of both the phase separation and the normal 
regimes. For sufficiently strong SOC, the phase separation region disappears, giving way to the 
superfluid state. We find that the tricritical point moves toward regime of low temperature, high 
magnetic field, and high polarization as the SOC increases. 
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Recent experimental realization of synthetic gauge 
field [H and spin-orbit coupling ^ opens a new arena to 
explore quantum many-body systems of ultracold atoms. 
The engineered SOC (with equal Rashba and Dresslhaus 
strength) in a neutral atomic BEC was accomplished by 
dressing two atomic spin states with a pair of lasers. The 
interaction between a quantum particles' spin and its mo- 
mentum is crucial for spin Hall effects [3] and topologi- 
cal insulators j^], which has captured great attention in 
the condensed matter community. The engineered SOC, 
equally applicable for bosons and fermions, allows for 
the realization of topological insulators and topologically 
nontrivial states [5|,|6[ in fermionic neutral atom systems, 
engendering broad interests in physical community. 

In anticipation of immediate experimental relevance 
involving SOC in fermionic atoms such as ^Li and 
intense theoretical attention has been paid to the 
physics of BEC-BCS crossover and polarized Fermi 



gases |lll - ll3l | in the presence of SOC . The SOC has been 



predicted to lead to various new phenomena. In particu- 
lar, for the two-body problem, it gives rise to a two-body 
bound state even on the BCS side (a^ < 0) of a res- 



onance 



14| . For the many-body physics at mean- field 



level, it enhances BCS pairing via the increased density 
of states at low energy, and leads to anisotropic super- 
fluids through mixing the spin singlet and triplet compo- 



nents 



Polarized fermionic condensates have been the focus 
of both theoretical and experimental research over the 
past years [31 . One of the key questions to ask is that 
how SOC reshape our understanding of these exciting 
systems. So far, most of the theoretical studies |8l-[l3| 
have focused mainly on zero temperature, leaving the 
physics at finite temperature which is experimentally rel- 
evant largely intact. In this paper, we are trying to ad- 
dressing this question by conducting the following stud- 
ies: Firstly, we map out the finite-temperature phase dia- 
gram at BCS side where mean-field theory gives quanti- 
tatively reasonable results. By determining the behavior 



of the tricritical point as a function of SOC strength, we 
can completely characterize the topology of the phase di- 
agram without recourse to an extensive numerical treat- 
Secondly, we examine the physics at uni- 
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ment 

tary, which is experimentally relevant and theoretically 
interesting. Specifically , we consider the effect of SOC 
on the "spin susceptibility" and the critical tempera- 
ture. Finally, we investigate the fate of breached pairing 
states [l^ under SOC through the correlation functions. 

We consider a homogeneous two species polarized 
Fermi gases interacting via an attractive contact poten- 
tial with an isotropic in-plane Rashba spin-orbit cou- 
pling, described by the following Hamiltonian, 

fd'rJ2 ^i(r)(^-M.)^.(r) 
/ ^-^ 2m 

-5 / Sxi,\(x)i,\(v)i,^{v)i,^(v) + Hso- (1) 



Here, Hso — A^j^ fcj^ [e^'''"'CktCk4. + /j-c] , with the 
transverse momentum k_L = (k^jky) and tp^ = Arg{kx + 
iky). The strength of spin-orbit coupling A can be tuned 
by atom-laser interaction Q . We define the chemical po- 
tential /i and "Zeeman" field h such that = fj, + h and 
fi^ — II ~ h. The spin imbalance between the two species 
is denoted by the polarization P = (n-f- — nj^)/(n-|- -|- nj,). 
We consider pairing between different hyperfine species 
of the same atom, so we restrict ourself to a single mass 
m. The interaction strength g is expressed in terms of 
the s-wave scattering length using the prescription: 
= + vEk2k' ■^here V is the volume and 
Ek = k^/2m (for convenience, we set h — ks ~ 1). We 
also define the Fermi momentum using kp = (Stt^tt,)-'^/'^, 
with total density n = n^ + ni, so that the Fermi velocity 
is vp ^ kp/m. Throughout our calculation, we will keep 
n fixed. 

Within the framework of imaginary-time field inte- 
gral, the partition function of the system is Z = 
/ c?[V', V'] exp (— S'[-0, ?/;]) with the action S'[-0,?A] = 



2 



/ [X^cr ^a^rfpa + H{tp, -0)] ■ Introducing a bosonic 
field A(r, t), which is believed to encapsulate the relevant 
low-energy degrees of freedom, we perform a Hubbard- 
Stratonovich transformation, then the action becomes 

/ dTcPr {^-^ - - AV'^V't) ■ To bring the 

action in a compact form, we define a four-dimensional 
vector ^fjc = (■(/'i(.-|.V'k;4,V'-ktV'-ki) ■ Then the action can 
be casted as S' = J dTj2k [^*k(-^^"^)*k + Ck] , with 
the inversed Green's function defined as 



V 



0_ 
A 



-A 






-A 

-dr^ Ckt 



A 





-dr + Cki/ 



with Ckt = Ck - + Cki = ek - /i - and Ck = Ck - A* 
. Integrating out the fermionic degrees of freedom, we 
obtain the effective action 

5e//= / dTd3r^-irrln(-g-i) + /3^ek. (2) 
J g z ^ 

Setting A = Ao -f (5A, and Q^"^ — t/^^|A=Aoj we can 
write = Q^"^ + S. Expanding the effective action to 
the second order in the fiuctuation S, we approximate 



the effective action to be Seff 



PYY^ SA{-q)SA{q) 



So 

- \Tr[go{k)j:{~q)go{k 



Sq, with 5*0 = 



, 5*3 — 

. 4 ..-"v^ mi)] ^ 

PVj: T-\q)A{-q)6A{q), where k = (k,*u;„), q = 
(q, iun) , and the quasiparticle excitation spectrum i?k± 
determined by El_j_ = + + h"^ + X^kj^ ± 

While i?k+ is always gapped. 



-Bk- accommodates gapless excitations distributed sym- 
metrically along fcz = axis at fc^ = 0: (1) for fj, ^ 
y^h'^ — Aq, it has four gapless excitation points at kz = 



±Jn±y/h^ - A2; (2) for ^ < y/h^ - A2, it only has 



A2 

^0- 



two gapless excitation points a.t kz = ^ ~ \Jh? 

Peculiar properties of the excitation spectrum are illus- 
trated in Fig. [1] where the isoenergy surface for iJk± — 
O.SEp at unitary (l/kpa = 0) is shown at zero temper- 
ature. The red dash curve is for i?k-, the blue solid 
curve is for -Bk-i-, and the green dash dotted circle is for 
a spherical Fermi surface. The isoenergy surface is sym- 
metric with respect to fc^ = 0, and possesses rotation 
symmetry along z-axis. For balanced superfluid (/i = 0), 

i?k± — \J (l^kl ± Afc^)^ -I- Aq, as depicted in panel (a) 
and (b), with A — 0.125?;^ and A = 0.25w_f respectively. 
The anisotropy of the isoenergy surface increases as one 
increases the strength of SOC A. It is interesting to no- 
tice that there exists two branches of isoenergy for both 
i?k+ and £'k- , due to the positiveness of chemical poten- 
tial in this case. This will lead to the enhancement of 
BCS pairing through increasing density of states around 
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FIG. 1. (Color online) Isoenergy surface (iJk± = O.S-Ef) 
for the quasiparticle excitation spectrum at unitarity where 
l/kpa = at T = 0: (a) = 0, A = 0.l2^vp- (b) /i = 0, 
A = 0.2^vf; (c) h = O.lEp; A = 0.1257;f; (d) h = O.lEp; 
A = 0.25wf. The red dash line is plotted for -Ek+, the blue 
solid line is for i?k-, and the green dash dotted circle is for a 
spherical fermi surface, plotted for comparison. 



Fermi surface. i?k+ and -Ek- merges at fcj^ = at which 
effects of SOC vanishes. For h — O.lEp, the isoenergy 
surface is shown in panel (c) and (d), with A = 0.125vp 
and A — 0.25vf respectively. Here we only have one 
branch for both £'k+ and i?k-- Interestingly, the curve 
for £'k+ develops a lunar structure for low A, as shown 
in panel (c), and increasing A makes it blunt. 

Phase diagram at finite temperature plays a key role 
in characterizing polarized fermionic condensates [l7l . [20| . 
For a fixed interaction strength l/kpa and a fixed SOC 
strength A, the phase diagram could be determined by 
the plane spanned by the temperature T = 1//3 and the 
Zeeman field h. At sufficient low polarization or chemi- 
cal potential difference h we expect a finite-temperature 
phase transition at which the superfluid order parame- 
ter vanishes continuously. Conversely, at low tempera- 
ture, superfluidity is destroyed in a flrst order fashion 
with increasing h. Across this phase transition at fixed 
h, the polarization jumps discontinuously. To determine 
the position of the phase boundaries, we must minimize 
the mean-field grand potential CIq — Sq / (3 with respect to 
the BCS order parameter Ao. Such a mean-field analysis 
should provides a qualitative reasonable description at 
weak-coupling BCS regime. In the absence of SOC (A = 
0), it is well known that there exists a finite-temperature 
tricritical point in the BCS limit, which is a natural con- 
sequence of having a first order transition from the super- 
fluid phase (SF) to the normal phase (N) at T = and a 
second order transition at zero polarization. First inves- 
tigated by Sarma f2l'| in the context of superconductivity 
in the presence of a magnetic field h, the BCS tricritical 
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FIG. 2. (Color online) Upper panel: finite-temperature phase 
diagram as a function of T and h at l/kpcis = —1 (BCS 
side). There are four different phases: the N state, the PS 
state, the SF state, and the magnetized superfluid (SFa/). 
Above the tricritical point, the transition line separating the 
broken-symmetry state (SFm) and the symmetric state (N) 
is of second order. Below the tricritical point (TP) it changes 
to the first order. Lower panel: the evolution of the tricritical 
point {Ttri/Tp, htri/Ep) as a function of SOC strength A. 

is located at (Tc^t / Aq, /icrit/Ao) = (0.3188,0.6061) 
where Aq = S/e^E'i? exp (— 7r/2|fcj?a|). The phase dia- 
gram spanned by T and h at 1/kpa = —1 for various 
SOC strength is shown in the upper panel of Fig. [21 It 
consists of four different phases: the superfluid state with 
zero polarization (SF), the magnetized superfluid state 
(SF7\/), the normal state (N), and the phase separation 
(PS) regime enclosed by the the first order line and the 
second order line. As the strength of SOC A increases, the 
area of phase separation region diminishes, and eventu- 
ally disappears for sufficient large A. With the increasing 
of SOC, the regime of SF diminishes very sharply, giving 
way to SF^f . The intersection of the second order line 
and the first order line gives the position of the tricritical 
point, denoted as TP in Fig. O The evolution of the tri- 
critical point is shown in the down panel. As A increases, 
Ttri decreases, while htri increases. 

The finite temperature phase diagram for l/kpag = 
— 1 spanned by T and P is shown in Fig. [3] The ef- 
fect of SOC is dramatic at zero temperature: the sys- 
tem is unstable to phase separation at any polarization 
without SOC; however, as the SOC is turned on, the 
system is a "magnetized" superfluid, in which the super- 
fluid component and the normal component coexist in 
an isotropic and homogeneous fashion. At finite temper- 
ature, by increasing the SOC strength A, both the regions 
of normal and the phase separation diminishes, leaving 
the broadened regime of SF. The tricritical point moves 
towards high polarization and low temperature when A 
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FIG. 3. (Color online) Finite-temperature phase diagram in 
the plane of T and P at l/kpaa = —1. The inset shows the 
corresponding polarization Ptri for the tricritical point as a 
function of SOC strength A. The phase SF is along the line 
of P = 0. The notation is the same as in Fig. [S] 



is increased, as seen in the inset. 

The effect of SOC on the stability of the system man- 
ifests itself on the spin susceptibility of the system. We 
examine the stability of the superfluid phase at unitarity 
at zero temperature where mean-field theory should give 
qualitative reasonable arguments. As shown in Fig. |31 at 
A = 0, the system is unstable to phase separation as the 
slope of SP/Sh is always either zero or negative. There 
exists a critical polarization (here Pc = 0.68), above 
which the system reverts to the normal state. When 
the SOC is turned on, the superfluid state could support 
both low polarization and high polarization, in contrast 
to what we saw Fig. [2] at 1/kpa ~ —1 where the super- 
fluid state only support low polarization. At a sufficient 
large A, the slope of the whole curve becomes positive, 
indicating that it is able to sustain any polarization. 

One of the most important questions to ask about 
the unitary superfluid is that how the critical temper- 
ature varies with SOC strength. At flnite temperature, 
the contribution from non-condensed pairs to the density 
n = dfl/dfi becomes important. This contribution is nec- 
essary to approach the transition temperature of an ideal 
Bose gas in the molecular limit where Tbec = 0.218Tp 
and can be included in the non-condensed phase (Aq = 0) 
through the gaussian contribution to the grand potential: 
fig = (l/2/3T^) I]q^.,j,^ lnr"i(q, w™), where Vra denotes 
the bosonic Matsubara frequencies. For numerical con- 
venience, we treat the fluctuation by adopting the NSR 
scheme [23| . The critical temperature for a balanced su- 
perfluid is shown in Fig. |4] at unitarity. The critical tem- 
perature calculated from mean-field theory Tco, starts at 
a high value and increases slowly at small SOC. Taken 
account of the gaussian fluctuation, the critical tempera- 
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FIG. 4. (Color online) Left panel: the polarization P — 
I as a function of maenetic field h for various SOC 
strength A at zero temperature at unitarity. Right panel: the 
critical temperature for balanced superfluid at unitarity: Tco 
is calculated from mean-field theory and Teg is calculated by 
taking account of the Nozieres Schmitt-Rind correction. 



ture starts at about 0.224Tf and increases almost linearly 
with increasing SOC strength A for low A, hinting at a 
possible way of realizing high- Tc superfluids. 

Another interesting question concerning two-species 
spin imbalanced Fer mi g ases is how the picture of 
breached pairing state jl9l| get modified. Breached pair- 
ing is characterized by a phase separation in momen- 
tum space between the excess of majority species t and 
the minority species J, in the superfluid state. Signa- 
tures of phase separation are visible in the momentum 
distribution, nko-, and correlation function C^^(k) = 
I < V'-kj.V'kt > 1- Requiring the pairing amplitude to 
be zero, one finds the region for phase separation: (1) 

fci = 0; and (2) \k^\ e [0, + + Al] if ^ 

y^h^ + A^ or \k, \ e [y^M- V^^^+Af, ^Jfi+^h^T^ol 
if ^ > ^/hF+'^. Referring to Fig.[5l for P = 0.7 shown 
in panel (a) and (c) , there exists two typical momenta kd 
and kc2 between which the minority species J, is depleted 
and the majority species has full occupation, reminiscent 
of breached pairing state with two Fermi surfaces (i3P2); 
while for P — 0.9 shown in panel (b) and (d), there 
exists a typical momentum kc below which the minority 
species is depleted and the majority species becomes fully 
occupied, reminiscent of BPi. In both cases, the correla- 
tion function shows a "hole" for momenta less than the 
Fermi momentum of the majority quasiparticles. The 
momentum distribution and the pairing amplitude bear 
consequences for experimental observation. The single- 
particle momentum distribution of trapped Fermi gases 
is routinely observed by timc-of- flight measurements ,24] . 

To summarize, we have identified a series of new fea- 
tures arising from spin-orbit coupling. We hope that cur- 
rent work will add new excitement to the surging field 
of cold atom physics involving artificial gauge field and 
spin-orbit coupling. 
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FIG. 5. (Color online) The momentum distribution Tiko- and 
correlation function C^|(k) — \ < -i/j-kj.V'kt > I a-t unitarity 
at zero temperature with SO coupling strength A = 0.2vf for 
two typical polarization: P = 0.7 (left panel); P = 0.9 (right 
panel) . 



We are grateful to Han Pu, Qi Zhou and Fei Zhou for 
helpful discussions. This work was supported by NSFC 
under grants Nos. 10934010, 60978019, the NKBRSFC 
under grants Nos. 2009CB930701, 2010CB922904, 
2011CB921502, 2Q12CB821300, and NSFC-RGC under 
grants Nos. 11061160490 and 1386-N-HKU748/10. 



[1] 
[2] 

[3: 
[4: 
[s: 
[6: 

[7 

[8 

[9 
[10 

[11 

[12 
[13 
[14 



rliao08@gmail. com | 
yyx@iphy.ac.cn 

Y.-J. Lin, R. L. Compton, K. Jimenez- Garcia, J. V. 
Porto, and I. B. Spielman, Nature 462, 628 (2009). 
Y.-J. Lin, K. Jimenez-Garcia, and I. B. Spielman, Nature 
471, 83 (2011). 

D. Xiao, C. M.-C, and Q. Niu, Rev. Mod. Phys. 82, 1959 
(2010). 

M. Z. Hasan and C. L. Kane, Rev. Mod. Phys. 82, 3045 

(2010) . 

M. Sato, Y. Takahashi, and S. Fujimoto, Phys. Rev. Lett. 
103, 020401 (2009). 

J.D. Sau, R. Sensarma, S. Powell, LB. Spielman, and 

S. DasSarma, Phys. Rev. B 83, 140510 (2011). 

J. P. Vyasanankere, S. Zhang, and V. B. Shenoy, 

arxiv:1104.5633. 

M. Gong, S. Tewari, and C. Zhang, arxiv:1105.1796, to 
appear in PRL. 

Z.-Q. Yu and H. Zhai, arxiv:1105.2250. 

H. Hu, L. Jiang, X.-J. Liu, and H. Pu, arxiv:1105.2488, 

to appear in PRL. 

M. Iskin and A. L. Subasi, Phys. Rev. Lett. 107, 050402 

(2011) . 

L. Han and C. A. R. S. Melo, arxiv:1106.3613. 

W. Yi and G.-C. Guo, Phys. Rev. A 84, 031608 (2011). 

J. P. Vyasanakere and V. B. Shenoy, Phys. Rev. B 83, 



5 



094515 (2011). 
[15] L. P. Gor'kov and E. I. Rashbar 87, 037004 (2001). 
[16] S. Giorgini, L. P. Pitacvskii, and S. Stringari, Rev. Mod. 

Phys. 80, 1215 (2008). 
[17] M. M. Parish, F. M. Marchetti, A. Lamacraft, and B. D. 

Simons, Nat. Phys. 3, 124 (2007). 
[18] M. M. Parish, F. M. Marchetti, A. Lamacraft, and B. D. 

Simons, Phys. Rev. Lett. 98, 160402 (2007). 
[19] W. V. Liu and F. Wilczek, Phys. Rev. Lett. 90, 047002 

(2003). 



[20] Y. il Shin, C. H. Schunck, A. Schirotzek, and W. Ketterle, 

Nature 451, 689 (2008). 
[21] G. Sarnia, J. Phys. Chem. Sohds 24, 1029 (1963). 
[22] R. Casalbuoni and G. NarduUi, Rev. Mod. Phys. 76, 263 

(2003). 

[23] P. Noziefes and S. Schmitt-Rind, J. of Low Temp. Phys. 

59, 195 (1985). 

[24] C. A. Regal, M. Greiner, S. Giorgini, M. Holland, and 
D. S. Jin, Phys. Rev. Lett. 95, 250404 (2005). 



